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Abstract 



We analyze K — > ttttj decays in the framework of Chiral Perturbation Theory. We study 
the different Dalitz plot distributions, trying to find regions where o(p 6 ) contributions 
could be more easily detected. To fulfill this program we compute all the the o(p 4 ) loop and 
counterterm contributions, finding a substantial agreement with the existing calculations 
and adding some small missing terms in K$ — > 7r + 7r~7. 
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1. Introduction 



The amplitudes for K — ► ttttj decays can be generally decomposed as the sum of two terms: the 
internal brcmsstrahlung (Ajb) and the direct emission (Abe) [1-4]. The first one, which represents only the 
bremsstrahlung contribution of the external charged particles, is completely predicted by QED in terms of 
the K — > tttt amplitude [5] . The second one, which is obtained by subtracting Ajb from the total amplitude, 
depends on the direct K — > irirj coupling and furnishes a test for mesonic interaction models. 

The pole for the photon energy going to zero would tend to enhance Ai b compared to Ad e ■ Nevertheless 
in Kl — > 7r + 7r~7 and K — > 7r ± 7r°7 decays the inner bremsstrahlung is respectively suppressed by the CP 
conservation and the AI =1/2 rule. As a result in these channels it is easier to study the direct emission, 
testing low energy interaction models. 

Chiral Perturbation Theory (ChPT) is supposed to be a reliable framework to study K — > ttttj decays 
[6, 7]. In this effective quantum field theory, based on symmetry principles, the physical amplitudes are 
expanded in powers of meson masses and momenta. At the lowest non-trivial order (o(p 2 j), only the inner 
brcmsstrahlung appears, while the direct emission, which requires more derivatives, starts at o(p A ). Experi- 
ments seem to indicate that in Kl — » 7r + 7r~7, and possibly in K ± — ► 7r ± 7r°7, not only o(p 4 ) but also o(p 6 ) 
magnetic contributions might be important [8, 4]. More in general, we can say that the question of the size 
of o(p 6 ) contribution in radiative non-leptonic kaon decays is very interesting, as it as been emphasized in 
the case of K L — > 7r°77 [9, 10]. 

In this paper we have analyzed K — ► nirj decays, looking for kinematical regions where order p 6 
electric contributions could be relevant and possibly detected. To this aim we have computed all the one 
loop amplitudes and the corresponding o(p i ) counterterms contributions. We substantially agree with the 
existing calculation for Kl — > 7r + 7r~7, — > 7r ± 7r°7 [4] and Ks — > ir + n 7 [11]. In addition we have 
calculated the small missing terms due irK and K r\ loops in Ks — > 7r + 7r~7. Using these results we find some 
kinematical region, in particular for Kl — > 7r + 7r _ 7 and possibly for K ± — > 7r ± 7r°7 decays, where order p 6 
electric contributions could be relevant. 

The paper is organized as follows: in section 2 we give a short discussion on kinematics, Low theorem 
and ChPT. In section 3, 4, 5 and 6 we analyze respectively Kl — > 7r + 7r _ 7, — ► 7r + 7r~7, — > 7r ± 7r°7 
and K L .s 7r°7r°7 decays. Section 7 contains some concluding remarks and finally in the appendices we 
discuss the loop functions and some variable transformations. 



2. Kinematics, Low theorem and ChPT 



Due to Lorentz and gauge invariance we can define two invariant amplitudes for the processes K(p K ) — > 
771(^1)^2(^2)7(9, e): the electric and the magnetic one. Using the dimensionless amplitudes E and M, defined 
in Ref.[4], we can write: 

A(K -> 7T7T7) = e»{q) [E{ Zl ){p iq p 2 ^ -p 2 qpi^) + M{z l )e^ pa p'ip p 2 q a ] /m 3 K , (1) 

with 

Zi = — o", Z 3 — — j - , Z3 — Z\ + Z2- (Z) 

m K m K 

Summing over photon helicities there is no interference among electric and magnetic terms: 
<9 2 r m K 



d Zl dz 2 (4tt) 3 



(\(E( Zl )\ 2 + \M(z t )\ 2 ) [ Zl z 2 (l - 2z 3 -r\- r\) - r\z\ - r\z\\ . (3) 



As we have already said in the introduction, generally one prefers also to decompose the total amplitude 
as a sum of inner bremsstrahlung and direct emission: 

A(K -» 7ri7r 2 7) = Ajb + A DE - (4) 
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Ajb, which in the classic limit would correspond to the charged particle radiation, in QED is completely 
predicted by the Low theorem [5] , which relates radiative and non-radiative amplitudes in the limit of photon 
energy going to zero. For K — > nnj transitions the theorem reads: 



A IB (K - 7T7T7) = e ( — - — I A(K 
\qPb qpa 



7T7T 



(5) 



where (p a ,Pb) = {p+,P-) for the neutral kaon decays and {p a ,Pb) = (p+,Pk) or (p a ,Pb) = (px,P-) for the 
charged ones. Thus in the E(zi) amplitude the two contributions Ajb and Af}^ ctrlc can interfere, differently 
from the M{z{) amplitude, where only a direct emission contribution appears: 

\A(K -> ^7r 27 )| 2 = \A IB \ 2 + 2 • Ue{A* IB A E D % ctnc } + \A E D % ctnc \ 2 + \A^ gneUc \ 2 . (6) 

By opportune kinematic integration, eq. (3) can be transformed in a photon energy distribution (see appendix 
B) . In the limit of the photon energy (E* in the kaon rest frame) going to zero one can write: 



dT a n 



(7) 



The content of Low theorem is a prediction for a and (5, from pure QED, in terms of the physical amplitude 
A(K^itit). 

We will now proceed analyzing K — > ttitj decays in the framework of ChPT. In this effective quantum 
field theory the octet of the pseudoscalar fields^ , 



tt°/V2 + v/V6 
K~ 



7T+ K+ 
-TT°/V2 + T]/V6 K° 

K° -vVys 



(8) 



is identified with the octet of pseudo Goldstone bosons coming from the spontaneous symmetry breaking of 
G = SU(3)l x SU(3)r into H = SU(3)l+r- The lagrangian, expressed in terms of 0, is expanded in powers 
of meson masses and momenta, and is dictated only by the transformation law of the interactions. 
At the lowest order p 2 , the strong part of the lagrangian is completely fixed and is given by: 



,(2) 



f {tv{D^D»U) + tr(*tft + f/ x t)} 



(9) 



where U — e 1 ^ 2 ^^ is a unitary 3x3 matrix which transforms linearly under G and \ is proportional to a 
scalar field, transforming like U , which must acquire a non- vanishing expectation value in order to reproduce 
the correct values of the meson masses [7]. At this order we can assume Fo = F n = Fk — 93MeV. The 
covariant derivative D^U is given by: 



(10) 



where Z M and r M are the external gauge fields of SU(3)l x SU(3)r. The electromagnetic field can be 
introduced setting \ jl = r jl = eQA^, where Q = diag(2/3, —1/3,-1/3). 

The |A5| = 1 weak lagrangian transforms under G as an (8l, 1_r) or a (27l, 1r)- At the lowest order 
p 2 can be written as: 



Cl 2 ^^ - G 8 F*tv (XD^D^U) 



-|AS| = 



+ G 21 Fl 



(U^D^U) 23 (D^U)n + -(U^D^U) 21 (D^U) 13 



h.c, 



(11) 



^ The chosen phase convention does not satisfy the Condon-Shortley-De Swart one [12]. 
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where A = (Xq — i\ 7 )/2 and A6,7 are the usual Gell-Mann matrices. Experimentally, from K — > tttt decays, 
the coupling constant G$ and G27 are fixed to be: 

|G 8 |^9-10- 6 Gey- 2 , f^-A" (12) 

At order p 2 , in K — > 77777 decays, since there are not enough powers of momenta, only an inner bremsstrahlung 
amplitude will appear. In agreement with the Low theorem we found: 

A IB (K -> 77777)°^) = e ( ^ - ^ - tttt) ^ 2 ). (13) 

Actually we will take eq.(13) as a definition of the inner bremsstrahlung amplitudes, meaning that this 
relation must hold order by order in ChPT. We recall also the o(p 2 ) results for the CP conserving K — > 7777 
amplitudes: 

A(K S - tt+tt-) ^ 2 ) = + 2(G 8 + \c^)F v {m\ - m 2 ), 

r (14) 

A{K+ 77^77^)°^ 2 ^ = + |G 27 f;« - m 2 ). 

ChPT is not a renormalizable theory but the requirement of unitarity implies that meson loops have to be 
considered. The loop contributions, which start at order p , introduce divergences which must be re-absorbed, 
order by order, by corresponding counterterms. While the physical amplitudes (loops + counterterms) are 
scale independent, the coefficients of the counterterms depend on the renormalization scale of the loops. 
Gasser and Leutwylwer [7] have classified the p operators for the strong and electromagnetic lagrangian 
and have determined their coefficients by comparison with the experiments. The authors of Ref.[13] have 
shown that these coefficients are well reproduced by vector and axial vector meson exchange choosing the 
renormalization scale \x around the p-mass. 

The effective Lagrangian at next to leading order, both in the strong and weak sector, has anomalous 
and non-anomalous operators [14]. At this order the anomalous operators contribute only to the magnetic 
amplitudes M(zi), which have been extensively studied in Ref. [3, 4] and will be not re-analyzed here. The 
non-anomalous octet part of the p 4 \AS\ = 1 Lagrangian, relevant for K — > 77777 decays, can be written as: 



4as|=i = G & F 2 [N 1A W 1A + N 15 W 15 + N W W 16 + N 17 W 17 ] + h.c, (15) 



where 



W14 


= itr 


w 15 


= itr 


w 16 


= itr 


w 17 


= itr 



(16) 



A = (A 6 - i\ 7 ) /2, u 2 = U, TV = iu^D^Uu^ = u],, /£" = uF^v) ± u^F^u (17) 



Ff = d^l v -d u F + i[l",r], Fg v = d»r v -0V + i[rV]. (18) 



and 

1 L - u 1 u 1 T ,i j, ± R 

In the electric transitions of K — > 77777 decays only the following counterterm combination appears: 

N m - N 15 - N w - N 17 . (19) 

Chiral symmetry alone does not predict the values of the constants iVj, these however can be estimated 
using various Vector Meson Dominance (VMD) models. Following Ref. [15], the combination in eq.(19) can 
be expressed in terms of the parameter kj: 



= (4tt) 2 [iV 14 - iV 15 - iV 16 - N 17 ] = -k f ^0- - o{l) (20) 
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At o{p & ) obviously the number of E-type counterterms increases, normalizing as and we 
can write in general: 

^ ) S |=i = G8 E JV i 6)w f )+h - c - (21) 

i 

We will not try to classify all the but we will consider only two typical operators: 

W e? = (^jl tr ( uAut (ux^uujru„ + u^mW)) , (22) 

< } = (4^4 tr ( uAut W> VpK, «,]]}) , (23) 

(Vjo is the "covariant" derivative defined in ref.[13], which for our purposes acts like d p ). The first one, 
similar to the one considered in Ref.[16], generates the same kinematical dependence of the o(p 4 ) operators, 
i.e. give a contribution to the electric amplitudes of K — > ttttj decays proportional to the one of the 
o(p ) counterterms. The second one, with more derivatives, generates a substantially different kinematical 
dependence. The size of the coefficients of these operator is an open question. With the normalization 
in eqs.(22)-(23), naive chiral power counting would suggest N^P <~ ~ o(l) nevertheless, for different 

arguments, sensibly larger estimates for and were made respectively in Ref.[16] and [17]. Waiting 
for further theoretical insight, we will not commit ourself to any particular model to predict the size of these 
coefficients, but we will investigate the possibility to bound or even to measure them. 



3. Kl — > 7T + 7T 7 



In this decay the electric and magnetic contribution of eq.(3) are competing and they have been both 
measured [8]. The inner bremsstrahlung violates CP through the Kl — > 7r + 7r~ amplitude and, neglecting 
the A J = 3/2 suppressed terms, can be written as: 

_ . . eeA(Ks — > ir + 7r~) , . . 

Eib (zi) = + r - L , Pi = P+, P2 = P- , 24 



which at the lowest order in ChPT is: 



f q(p 2 )^ n _ 2eeG 8 F 7r (m|- - m|) 

ib \ z i) ^TTr~7~\ — • ^ ' 

mK(z+z-) 



To decrease the theoretical uncertainties it is wiser to use in eq.(24) [11] the experimental value for \A(K S — ► 
7r + 7r _ )| [18] and the phases according the most recent analysis based on very general assumptions like 
analyticity, unitarity and crossing symmetries [19, 20]. Of course one has to use the phase convention of 
eq.(ll). 

The direct emission contribution can be generally decomposed in a multipole expansion [2]: 



*_) + ... 



E D e =Ei(z 3 ) + E 2 (z 3 ) ■ (z+ i ... 

M=M 1 (z 3 ) + M 2 {z 3 )-(z + -Z-) + ... ( ' 

We remind that even electric and odd magnetic multipoles are CP even, while the others are CP odd. 
The multipole expansion is particularly justified in these decays by the smallness of the phase space, which 
implies \z + — z_| < 0.2. Furthermore, local contribution to higher multipoles are suppressed by chiral power 
counting. So far only the inner bremsstrahlung and the dipole magnetic contribution have been observed. 
The results for the branching ratios are [8] : 

BR(K L -» 7r+Tr- 1 ) IB .E'>20 MeV = (1.32 ± .05)10- 5 , 

j. r 27 

BR(K L — > 7T 7T 7)l)£;,£;*>20 MeV — 

(2.95 ± .11)10"°. 
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For the IB both the spectrum and the rate are in agreement with eq.(24). The dipolc magnetic contribution 
can be written in first approximation as 



M 1 (z 3 )~M 1 -(l + cz 3 ), (28) 

where M" and c are fixed by the rate and the spectrum respectively. The large experimental vale for the slope: 
c exp = (— 1.7 ± 0.5) [21], which vanishes at o(p 4 ) in ChPT, carries an important dynamical informations. It 
is not only a pure o(p 6 ) contribution, but it is also substantially larger of a typical vector meson exchange 
effect [4], showing that, at least in the weak anomalous sector, local p 6 contributions are sizable and can 
play a role in discriminating different models. 

Regarding the electric multipoles, at o(p A ) there are no local CP conserving counterterms and conse- 
quently the loop contributions are finite [22, 23]. We have computed these amplitudes confirming the earlier 
results of Rcf.[4]. The potential contributions are given by irir, irK, Krj and KK intermediate states (see 
fig.l). Nevertheless, due to CP conservation, only the small ttK and Krj loops are non- vanishing. Their 
contributions can be parametrized using the functions h W K and hxr) defined in the appendix A: 

E\tlM) = eGsmK 8 { fg- ml) [hMz-) + h Kn (z-) Kk{z + ) - h Kv {z + )\ . (29) 

The functions hij depend only on one kinematical variable. Considering the Taylor expansion of these 
functions (hij(zi) = T^ n ^ij z i)^ w ^h tne coefficients reported table 1, we found: 



Ez^ ^t^:^ ^-- -- ,] - 



Consistently with the CP conservation only an E 2 term survives but it is suppressed by several factors 
compared to Eib in eq.(25): i) the angular momentum barrier, ii) the absence of the photon energy pole, 
iii) the chiral suppression factor ( 4 ™^ ) 2 and iv) the small coefficient 0.05 generated by the chiral loops. 

Other Ei contribution can arise at o(p 6 ). The operator is forbidden by CP conservation, like the 

o{p i ) counterterms, on the other hand the operator W E 6 ^ in cq.(23) gives: 

= +^<W^>. ( 31 ) 

Since <~ o(l) the E 2 multipole should be dominated by the local o(p 6 ) counterterms instead of the o{p i ) 
loop contributions. This is one of the main message of this paper. At higher order chiral loops will generate 
a small rescattering phase 5e 2 to the term in eq.(31). 

Now the question is whether the E 2 multipole can be detected. The electric contribution to eq.(3) is 

\Eib + E DE \ 2 ~ [\E IB \ 2 + 2Ke{E* IB E E2 )] 

= " " ^2 r 2 r 2 — t 1 + C int Z+Z-{Z+ - Z-)] , 

lib jy- 6j\ £/ 



where 

m 4 K cos((p e + 5$ -Se 2 ) at (6) 
48|e|7r 4 F4 

(5g ~ (39 ± 5)° [19] is the the (7r+7r-)/ =0 phase shift, £ ~ (43.67 ± 0.14)° [24] is the phase of e and Se 2 is 
expected to be small. Unfortunately the argument of the cosine in eq.(33) tends to be about 80°, smearing 
somewhat the effect. Nevertheless we can try to study the interference term C^ nt , and thus N E 6 2 ', measuring 
the (indirect) CP violating asymmetry under pion exchange, i.e. the observable: 

L N(K L -> Tr + n- T: cos(6») > 0) - N(K L -> 77+77-7; cos(6») < 0) 



N(K L -» 7r+7T^7; cos(0) > 0) + N(K L -> 77+77-7; cos(6>) < 0) ' 



(34) 



6 



where 9 is the angle between 7 and tt + momenta in the di-pion rest frame (see appendix B). For E* > 20 
MeV 

|Af B;>20 MeV) I - 1-2 • lCT 4 |Ct t | ~ 10- 3 |JVg|, (35) 

thus at least 10 11 /Cl are necessary in order to observe a non vanishing effect. 

We finally note that A L could receive a contribution also from the interference between E 2 and a direct 
CP violating multipolc E\. Nevertheless this contribution it is not enhanced by the interference with the 
inner bremsstrahlung and therefore should be suppressed respect to the previous one. 



4. Kg — > 7T+7T 7 



In this channel we will consider only the electric amplitudes since the dipole magnetic contribution is CP 
violating while higher multipolcs are phase space suppressed and do not interfere with electric amplitudes 
for unpolarized photons. The CP conserving bremsstrahlung term 



E IB {z l ) = + 



eA(K s — > 7r+7r~) 



m K (z + zJ) 

dominates giving a contribution to the branching ratio 



(Pi =P+, P2 =P-), 



BR(K S 



7T + 7r~ 



7) 



IB,E*>20 MeV 



4.80-10" 



(36) 



(37) 



and a Dalitz Plot distribution, reported for instance in Ref. [25, 2], which agree very nicely with the data [8]. 

At the lowest order in ChPT only the IB is present, the electric dipole component appear at the next 
order. The two gauge invariant contributions Eib and Ede are separated according to the ChPT definition 
of IB (eq.(13)). At order o(p 4 ) there are loops and counterterms. Since the counterterm contribution to 
A(K S — » tt + tt~^) D e 



4eG 8 m 3 K 



[N 14 - N 15 - N 16 - N 17 ] = 



eG s m 3 K (4) 



(38) 



is scale independent [22, 23], the direct emission loop contribution has to be finite [11]. As in the Kl case 
the possible intermediate states appearing in the meson loop are tttt, ttK, Kr\ and KK (see fig.l). In Ref. [11] 
only the tttt loop was computed, since this is the only one which has an absorptive part. Interestingly a large 
dispersive contribution was found, even larger than the absorptive one in the allowed kinematical regions. 
This of course urged a calculation of the other dispersive contributions and a confirmation of the previous 
7T7T loop calculation. Analogously to eq.(29), our result for the loop contribution can be written as: 



eG 8 m K (m 2 K - ml) 
87^ 



[4:h wv (-z 3 ) + K K {z+) + h Kv (z+) + (z + <-> z- 



(39) 



where h^n , k and Hk v are defined in the appendix A (the KK contribution vanishes in the limit m K + = 
rriKo). Regarding the irir loop we confirm the interesting result that the dispersive contribution is larger 



than the absorptive one, however we disagree on the analytic structure. We find indeed 



EW{E1) = - U< rnl){m% 2E*m K ) 



8ir 2 F v E* 2 



-E*m K (2ml - 2m 2 K ) + m\(m 2 K - m\) 



phi 



In 2 



(3-1 



In 2 



■A) In 



i + A) 



1 + /3 



(40) 



while the authors of Ref. [11] for the polynomial term had (3m 2 — 2m 2 K ). Though the difference is numerically 
insignificant has some deep meaning: i) the full amplitude is proportional to [rv? K — to 2 ) , i.e. to the weak 
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vertex of fig.l with the pions on-shcll, ii) the E^ op (E*) amplitude vanishes in the limit E* — > 0, as required 
by Low theorem'''. Of course with the help of these considerations one understands better why the KK loop 
contribution is zero. Furthermore it is also confirmed by our calculation that the irir loop contribution is 
dominant as suggested in Ref.[ll]. 

Another interesting point of our calculation is the kinematical dependence of the loop contributions. 
Using the Taylor expansion of the function hij in eq.(39) we obtain: 



EZ(zi) 



eG 8 m K (m 2 K - m|) 



[ 1.3 + l.Lz 3 - i(0.6 + 0.9z 3 ) 



(41) 



Since the kinematical dependence of eq.(41) is mild, in principle it could be possible to distinguish particular 
o(p 6 ) contributions, like the one of the operator W^: 



EfAzi) = - 



eG a m 3 K 
24tt 4 F3 e - 



(6) 



(42) 



On the other hand it is impossible to distinguish among the flat contributions of the o(p 4 ) counterterms and 
the operator W^: 



eG 8 m 5 K (6) 
327r 4 F 3 Si' 



(43) 



Of course all these contributions give a very small intcrfcrcncial branching ratio compared to the IB in 
eq.(37). As explained in the discussion regarding the bremsstrahlung of Kl — > tt + tt^j, for A(Ks — » ir + n~) 
we will use the experimental value of \A(Ks — » 7r + 7r~)| and Sq ~ (39 ± 5)° [20, 19]. Neglecting the p 6 
contribution, using the central value of <5° and varying we obtain the results reported in table 2. 

Differently from the K L case, due to the large IB, in K L — > tt + tt~j decays it is practically impossible 
to measure a CP violating interference between an even and an odd multipolc. 



5. K± 



7T ± 7r°7 



Since the initial and the final states of K ± — > 7r ± 7r°7 decays are not CP eigenstates, in these channels 
the inner bremsstrahlung amplitude together with the lowest electric and magnetic transitions are present, 
also in the limit of CP conservation. Actually Ejb is suppressed by the AI =1/2 rule: 



Ei B {zi) = 



±eA(K ± -» t^tt") 
m K {z±z s ) 



(Pi =P±, P2 =Po)- 



(44) 



Experimentalists generally choose one kinematical variable as the kinetic energy of the charged pion (T* in 
the CMS), which is not affected by the problem of identifying the ir° photons. In the range 55 MeV< T* < 
90 MeV the theoretical prediction for internal bremsstrahlung branching ratio is: 



BR(K ± — > 7T ± 7r°7)/B, 55MeV<T*<90MeV = 2.61 • 10" 



(45) 



After the subtraction of this contribution from the experimental branching ratio there is a clear evidence of 
a direct emission component [18]: 



BRiK* -> ^ ± ^°7)Z3B, 66MeV<T;<90MeV = (1.8 ± 0.4) • 10" 



(46) 



It is only one order of magnitude less than the IB in cq.(45), instead of the typical 10 2 <~ 10 3 suppression 
factor of the interferencial contribution, because both the electric and magnetic amplitudes are not suppressed 
by the AI = 1/2 rule. 



t The author of Ref.[ll] have recently checked their calculation confirming our result. 
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From the analysis of the Dalitz Plot distribution, the DE component seems more likely due to a magnetic 
transition [26], but the present data are not conclusive about this point. The theoretical discussions of the 
magnetic amplitudes has been done for instance in Ref.[3, 4]. Here we will concentrate on the electric ones 
in the framework of ChPT. At the lowest order, as usual, there is only the IB contribution: 



3m K (z±z 3 ) 



At next order there are loops and counterterms. The o(p 4 ) counterterm combination which appears is the 
same finite one of eq.(38): 

4 4 ^) = T^<\ (48) 

thus the loop contributions are finite. In the octet limit (G27 = 0) only the small nK and Kr\ loops are 
non- vanishing. We have computed them finding a little discrepancy with the earlier computation of Rcf . [4] . 
Our result is: 

( 4) , > eG 8 m K (m 2 K - ml) 
E ioop( z i) = T SZ2^ [Kk{z+) + h Kv (z+)} 



7 2, ( 49 ) 
eG 8 m K (m K - m v ) 

" ± 8^ [ " 16 + mZ+h 

weather in Ref.[4] the hxn function was multiplied by a factor 2/3*. Due to the smallness of the loop 
contributions this discrepancy is numerically irrelevant. 

As in the Kl — > ir + ir~j case, in — > 7r ± 7r°7 decays it turns out that the loop contributions are 
very small and thus the o(p 6 ) local operators may be relevant. However in this case the o(p 4 ) counterterms 
are not vanishing and should represent the dominant effect. To disentangle the o(p 6 ) contributions from 
the o(p 4 ) ones it is necessary to look for particular kinematical dependencies. As in the previous cases the 
operator gives only a flat contribution which renormalizes the one of the o(p 4 ) counterterms: 

On the other hand the operator gives a non-flat contribution: 



F (6), , eGsmx (6) 

^E 2 (*%) ~ T 247r 4 F 3 iV S 2 



3 



(51) 



which could be detected studying the T* dependence of the Dalitz Plot distribution (see appendix B). 



6. K L , S -» 7r u 7r u 7 

Since no charged particles are involved, $ — * 7r°7r°7 decays are completely different from those con- 
sidered before. The absence of charged particles implies the absence of the inner bremsstrahlung amplitude 
and thus a very small branching ratio. Furthermore these decays are suppressed by Bose statistics, which 
requires both electric and magnetic multipolcs to be even. As a consequence none of them has been observed 
yet. 

In the limit of CP conservation Kl — > 7r°7r°7 is a pure electric transition and Ks — > 7r°7r°7 is a pure 
magnetic one. We will consider only the former. As it has been shown in Ref.[27], in Kl — * 7r°7r°7 not 
only the o(p 4 ) local contributions are forbidden, but also the one-loop amplitude is vanishing. Therefore this 
decay can receive contributions at least of order p 6 in ChPT and could be very useful in order to fix, or even 



t The author of Ref.[4] have recently checked their calculation confirming our result. 
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to bound, these new couplings. Unfortunately it is very difficult to observe, due to the large background of 
K L -> 3tt°. 

A typical operator which gives a non vanishing contribution to Kl — > 7r°7r°7 is W^: 

Using eq.(52) we can parametrize the branching ratio of JTj, — ► 7r°7r°7 in terms of N^: 

i?i?(X L ^^ 7r° 7 )^10- 9 |iVf| 2 . (53) 

In Rcf.[17] the optimistic guess of BR(K L -> 7r°7r°7) ~ 1CT 8 , together with BR(K S -> 7r°7r°7) ~ 1CT 11 , 
was made. This prediction would suggest ~ 3, which is certainly bigger than the naive power counting 
estimate but cannot be excluded a priori. 



7. Conclusions 



The analysis presented here completes the previous o(p 4 ) calculation onif-> 7T7r 7 decays and explores 
the effects of o(p 6 ) local operators in the electric transitions. 

Regarding the o(p 4 ) DE loop contributions, we find substantial agreement with the existing calculations. 
In Ref.[ll] it was shown that the dispersive pion loop contribution is larger than the absorptive one. This has 
urged us to calculate also the ttK and r\K loops, which give only a dispersive contribution. We find that the 
naive expectation of Ref.[ll], that these contribution are small, is confirmed. Furthermore our calculation 
has clarified some erroneous considerations of Ref.[ll]. 

We have also analyzed the contributions to K — > irir-y decays of two typical o(p 6 ) operators (eqs.(22)- 
(23)), which have been studied for different purposes in the literature. Interestingly we find that in Kl — * 
7r + 7r~7 the local o(p 6 ) contributions can be substantially larger than the o(p A ) ones. 
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Appendix A 



Following Ref.[4], in order to write the loop functions hij(z) we first introduce the expansion of the 
three point integrals with q 2 = 0: 



dH 

[P - m 2 } [(I + q) 2 - m 2 ] [{I - p) 2 - m 2 



ig^C 2 o(p 2 Ap + q) 2 ,mlm 2 ) + V.g" terms". 



With this definition, the finite functions hij (z) are given by: 



(Al) 



hi^z) = (Att) 2 ^- [C w (p 2 ,{ P + q ) 2 ,mlrn 2 )-C w {p 2 ,p 2 ,mlm 2 )\ , 
pq 



(A.2) 



where z = pq/m 2 K and p 2 = m 2 (for h^K and hxr/) or p 2 = m 2 K (for h^^). 

These functions can be explicitly written in terms of one dimensional integrals. Defining: 



Hz) 

fc 



2 f 1 dx ( m 2 (1 - x) + xm 2 - x(l - x){p 2 + 2zm 2 K ) 



zm 



k Jo 



— In 

x 



2zm 2 K 



m 2 (l — x) + xm 2 — x(l — x)p 2 

I dx In 

Jo 



p 2 + m 2 — m 2 + 2zm 2 K f 1 i i f m 2 (1 - x) + xm 2 — x(l — x)(p 2 + 2zm 2 K ) 



(1 — x) + xm 2 — x(l — x)p 2 



2 r 1 



m 

V Jo 



dx In 



mf (1 — x) + xm 2 — x(l — x)p 



+ 



2 r i 



ml f 
P 2 Jo 



dx In 



m\ J p' 

mf (1 — x) + xm 2 — x(l — x)p 2 \ m 
m 2 J p< 



(A3) 



we can write: 



hv(z) = -{f 1 (z) + f 2 (z) + f c ]. 



(AA) 



As expected by the Low theorem, (fi(z) + f2(z)) ^> — f c , so that hij(z) has no pole for z — > . 
In the case m, = mj the integrals f\ , f 2 and f c can be done explicitly. For the rnr function (rtii = rrij 
m 7T, P 2 — m %) we found: 



K ^ z )=—{(l + 2z) 



(3\n 



1 + /3 



m 



K 



In' 



0-1 
1 + Po 

A>-i 



-A) In 



i + A) 



In" 



1 + /3 
(3-1 



2z}, 



(A5) 



where 



(3o 



Am 2 



and (3 = a 1 



Ami 



K 



m 2 K (l + 2z) 



(A.6) 



It is interesting to no note that the function h vv (z) is simply related to the function H(z), defined in 
Rcf. [28] , which appears in If — > 77* decays^ : 



h vw (z) = -H(l + 2z). 



(A.6) 



t We thank G. Ecker for clarifying us this point 



If 



Appendix B 



In Kl.s — * 7r + 7r~7 decays, the variables which can be more easily studied from the experimental point 
of view are: i) the photon energy in the kaon rest frame {E*), ii) the angle between 7 and tt + momenta in 
the di-pion rest frame (9). The relations between (E*, 9) and the Zi are: 



Z3 



E* 
m K 



^■(lT/3cos(0)), 



{B.I) 



where (5 = yjl — 4mJ/(m| — 2m kE*). The kincmatical limits on E* and 6 are given by: 

-1 < cos(6>) < 1. 



o<e; < 



Ami 



2m,K 



Finally the differential rate in terms of these variables is: 
d 2 T Elp d 2 Y 



dE*dcos(9) 2m 2 K dzidz 2 



\E\ Z + \B[ 



e; 3 /3 3 

5127r 3 r 



2_e; 

m K 



sin 2 (0). 



(B.2) 



(B.3) 



In K ± — > 7r ± 7r°7 decays the situation is different, it is more useful to study the differential rate as a 
function of [26]: i) the charged pion kinetic energy in the kaon rest frame (T*), ii) the adimensional variable 
W 2 = (qp K )(qp±)/(m 2 r+ m 2 ^). These variables are related to the Zi by: 



z± 



2m\ 



[m 2 K + ml+ - mlo - 2m K ra^ - 2m K T*] , z 3 z ± 



m 



2 vv ■ 



K 



While for T* it is easy to write the limits: 



< T* < 



(m K - m^+) 2 - mlo 
2mK 



(BA) 



(B.5) 



for W the expressions are cumbersome and we refer to the figure 3.3 of Ref. [2]. The advantage in using 
these variables lies in the fact that, through the VF 2 -dcpendcnce, it is easier to disentangle the different 
contributions of inner bremsstrahlung, direct emission and interference [29] : 



d 2 T 



d 2 T 



IB 



dT*dW 2 dT*dW 2 
where A = A(K ± ^ t^tt ). 



1 + 



l l±2Vte ( ^ W 2 



m K 



\1* 



+ 



"K 



Ede 



eA 



+ 



M DE 



eA 



(B.6) 
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Tables and Figures 



Table 1 

Coefficients of the Taylor expansion of the functions hij, for ij = nn, ttK and K r\. 



h° 

a ij 



HIT 

ttK 
Kr\ 



+ .36 - A5i 
-.12 
-.04 



.30 + .22i 
-.045 
-.005 



.37- A2i 

< io- 2 

< io- 2 



Table 2 

Inner bremsstrahlung and direct emission contributions to the branching ratios of K$ — > 7r + 7r~7 decay, 
for different values of the cut in the photon energy. Different values of the o(p ) counterterms are chosen. 



E* > 20MeV E* > 50MeV E* > lOOMeV 



10 3 


IB 




4.80 


1.75 


.31 


10 6 


DE (k f 


= 0) 


-6.2 


-5.0 


-2.0 


10 6 


DE (k f 


= +.5) 


-10.5 


-8.3 


-3.3 


10 6 


DE (k f 


= +1) 


-14.8 


-11.7 


-4.7 


10 6 


DE (k f 


--•5) 


-1.9 


-1.6 


-.6 


10 6 


DE (k f 


= -1) 


+2.4 


+1.8 


+.7 



Figure 1 

o(p 4 ) loop diagrams for K —> ttttj direct emission amplitudes. The symbols o and • represent the strong 
and the weak vertex, respectively. The photon has to be attached to any charged line or to any vertex. 
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